We compute the afterglow of gamma-ray bursts pro- 
duced by purely electromagnetic outflows to see if it 
shows characteristic signatures differing from those ob- 
tained with the standard internal/external shock model. 
Using a simple approach for the injection of electromag- 
netic energy to the forward shock we obtain the afterglow 
evolution both during the period of activity of the central 
source and after. Our method equally applies to a vari- 
able source. Afterglow light curves in the visible and X- 
ray bands are computed both for a uniform medium and a 
stellar wind environment. They are brighter at early times 
than afterglows obtained with the internal/external shock 
model but relying only on these differences to discriminate 
- - - between models is not sufficient. 
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1. Introduction 

Lyutikov and Blandford (2003) proposed an alternative to 
the standard fireball model where the central engine pro- 
duces a purely electromagnetic outflow instead of a rela- 
tivistic baryonic wind. Observationally this electromagne- 
tic model (hereafter EMM) differs from the standard in- 
ternal/external shock model by the absence of any reverse 
shock contribution, a different early afterglow evolution 
and a high polarization of the prompt emission (Lyutikov, 
2004). In this paper we concentrate on the early afterglow 
(while the central source is still active) and compare the 
EMM to the standard model in X-rays and the visible for a 
uniform external medium or a wind environment. In Sect. 2 
we obtain simple equations that govern the evolution of 
the forward shock propagating in the burst environment. 
Their solutions are used in Sect. 3 to compute afterglow 
light curves which are compared to those obtained in the 
standard model for the same total injected energy. We 
discuss our results in Sect. 4 and conclude that it likely 
will be difficult to decide between models from afterglow 
observations only. 



2. Dynamics of the forward shock 

In the context of the EMM, electromagnetic energy re- 
leased by the central source directly leads to the formation 
of a forward shock propagating in the external medium. 
We obtain the evolution of this forward shock by writing 
the conservation of energy-momentum of the swept-up 
mass as it accumulates electromagnetic energy 



Eem 
Eem 



Mc 2 = MTT lC 2 
: Mf3TT iC 2 



(1) 



where M is the swept-up mass, (3 and T the velocity and 
Lorcntz factor for the bulk motion of the shocked material 



and Ti the Lorentz factor for internal motions. The energy 
Eem received from the source reads 



Eem = 



L EM (t)dt 



(2) 



where LemW is the source electromagnetic power and t s 
the time in the source frame (which is also the observer 
time modulo the (1 + z) factor for time dilation). 

This can be easily understood if we define the time 
t e when the shocked material is located at radius R. The 
emitted afterglow signal will reach the observer at 



^-obs 



te 



(D-R) 



(3) 



where D is the distance of the source to the observer. At 
time t ei all power produced by the source up to a time t s 
such as 



c(t e -t s ) = R 



(4) 



has been received by the moving shell and Eq.(3) and (4) 
then lead to 



D 

ts + — - 

C 



tabs = t s H . (5) 

Going back to Eq.(l) and after elimination of Ti we get 

2Mc 2 r 2 ~ Eem (6) 

this approximation being valid in the relativistic phase 
(r 1). Differentiation of Eq.(6) with respect to observer 
time (with the simplified notation t = t Q b s ) yields 



£em(*) 



IQitAc 2 



3-s 

where we have written 

M(R) = i^-i? 3 -* 
3-s 



i?3- sr ^E + (3 _ s)ci? ; 



dt 



(7) 



(8) 
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Fig. 1. Evolution of the Lorentz factor in the different models considered; left panel: uniform external medium of 
density n = 1 cm -3 ; right panel: stellar wind with A* = A/(5 10 11 g.cm -1 ) = 0.1. The full, dotted and dashed lines 
respectively correspond to the EMM with a constant Lem = 10 52 erg.s -1 , a variable Lem(£) given by Eq.(13) and to 
the standard internal/external shock model with E = 10 52 erg.s -1 and an initial distribution of the Lorentz factor in 
the relativistic wind given by Eq.(14). 



with A = p and s = for a uniform medium of density p 
and ^4 = 5 10 n ^4* g.cm -1 and s — 2 for a wind environ- 
ment. With the additional relation between observer time 
and shock radius 



dR 
dt 



= 2cT 2 



(9) 



the problem can be solved for any law Lem(*)- With a 
constant Lem it can be easily shown that the solutions of 
Eq.(7) and (9) are 



r = (Q/2) 1/2 (ci)- 1 / 4 and R = 2Q{ct) 



1/2 



with 



3L 



EM 



1/4 



V 32wpc 3 , 
for a uniform medium and 

1/4 
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and R = 2cT A t 



(11) 



(12) 



for a stellar wind. With a variable £em(*)) Eq.(7) and (9) 
are integrated using Eq.(10) or (12) as initial conditions 
with L E m = L EM (t = 0). Figure 1 shows the resulting 
T(t) for a constant L E m 

by 



10 52 erg.s 1 or L EM (t) given 



£bm(*) = 10 



">2 



1 + C0S27T 



10 s 



erg.s 



(13) 



In both cases the source is supposed to be active for 50 s. 
For comparison we also plot in Fig.l the evolution of the 
Lorentz factor computed in the standard internal/external 
shock model. In that case, we assume that a constant ki- 
netic power E = 10 52 erg.s -1 is injected in the relativistic 
wind for a duration i\y — 50 s. We consider a single pulse 
burst produced by a distribution of the Lorentz factor of 
the form 



(io) r(t) 



■ cos 



fok;) (14) 



if t < 0.2 t w and T(t) = T max if t > 0.2t w ; T max = 200 
and r m ; n = 50 are the maximum and minimum values 
of the Lorentz factor so that the most rapid part of the 
wind is decelerated by the slower one which was emitted 
previously. To follow the wind evolution and its interac- 
tion with the external medium we represent the flow with 
a large number (> 1000) of discrete shells which interact 
by direct collision only, pressure waves being neglected 
(Daigne and Mochkovitch, 1998). The action of the exter- 
nal medium is included by progressively adding the swept- 
up mass to the front shell of the wind. In all cases, the 
Lorentz factor relaxes to the standard Blandford-McKee 
(1976) solution after a few 10 2 s. However at early times 
r is larger in the EMM, especially for a uniform external 
medium. 
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Fig. 2. The two frequencies v m (full lines) and v c (dashed lines) in the EMM (thick lines) and the standard model 
(thin lines); left panel: uniform external medium of density n = 1 cm~ 3 ; right panel: stellar wind with A* =0.1. 



3. Afterglow calculation 

3.1. Method 

Using the dynamical evolution of the forward shock ob- 
tained above we have calculated afterglow lightcurves in 
X-rays and in the V band for both the electromagnetic and 
the standard model. We obtain the two critical frequen- 
cies v m and v c following Sari, Piran and Narayan (1998). 
They are represented in Fig. 2 for a uniform medium and 
a stellar wind environment. In the EMM with a uniform 
medium and a constant source, the ratio v m jv c is given 

by 

— = 0.065 £e._ie|._ 3 n L 52 (15) 

where n is the density in cm~ 3 , L 52 the electromagne- 
tic power in units of 10 52 erg.s -1 , e e ,-i = Ee/ICT 1 and 
e.B-3, = £b/10~ 3 (the numerical factor in Eq.(15) cor- 
responding to the electron power law index p — 2.5). 
The cooling regime therefore remains fixed as long as the 
source is active (and does not vary). In Fig. 2 we have 
adopted e e ,-i = £5,-3 = £52 = 1 and n = 1 or A* = 0.1 
in the uniform medium or wind case respectively. 

With these parameters the afterglow is always in the 
slow cooling regime in the uniform medium case while in 
the wind case it is in fast cooling while the source is active 
and moves to slow cooling shortly after. 

3.2. Results 

We have represented in Fig. 3 afterglow light curves in X- 
rays (2 - 10 keV) and in the V band for a uniform exter- 



nal medium and a stellar wind. The burst parameters are 
identical to those used in Fig. 2. For comparison we also 
show in Fig. 3 standard afterglows computed in the in- 
ternal/external shock model with the same total injected 
energy. A redshift z = 2, typical of SWIFT bursts has 
been adopted (Jakobsson et al., 2006). In all cases the af- 
terglows are brighter at early times in the EMM while at 
late times (in the Blandford-McKee regime) the models 
coincide. The difference in early evolution is larger for a 
uniform external medium, especially in X-rays. This is due 
to the large initial Lorentz factor in the EMM compared 
to the standard model (see Fig.l) which leads to a higher 
electron Lorentz factor in the shocked external medium. 
The frequency v m is then in the X-ray band already at 
very early times while it is still in the UV/visible in the 
standard model. In the wind case, the differences between 
the EMM and the standard case are smaller and the two 
models will be therefore more difficult to distinguish. 

For the EMM we have also calculated afterglows when 
the source is variable (with Lem(0 given by Eq.(13)). The 
results are shown in Fig. 4 for a uniform external medium 
and the first 200 seconds of evolution (the source being 
active for 150 s in the observer frame). The left panel 
in Fig. 4 shows the simple calculation which only includes 
line of sight emission. It produces a highly variable X-ray 
light curve but off-axis effects (time delays and spectral 
softening) will smear out any variability occurring on a 
time scale shorter than 



where R and T are the radius and Lorentz factor of the 
emitting shell. The right panel of Fig. 4 shows At as a func- 
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Fig. 3. Afterglow comparison between the EMM (solid line) and the standard model (dotted line). Top: X-ray band 
(2 - 10 keV); Bottom: V band; Left column: uniform medium of density n = 1 cm~ 3 ; right column: stellar wind with 
A* = 0.1. A rcdshift z — 2 has been assumed. 



tion of observer time together with the X-ray light curve 
now obtained with a detailed calculation including off-axis 
effects (Granot, Piran & Sari, 1999; Woods & Loeb, 1999). 
Except for the first two or three peaks which partially 
subsist because At is initially small, the rest of the curve 
becomes nearly monotonic. In the visible the variability is 
barely visible, even without including off-axis effects. 



4. Discussion and conclusion 

The lightcurves in Fig. 3 show that the EMM and the stan- 
dard model notably differ at early times (during the period 
of source activity). However relying on these differences 
alone to identify the physical origin of GRBs will be a dif- 
ficult task requiring a very early follow-up of the afterglow. 
In X-rays, SWIFT should be able to do that (at least in 
some cases) but the problem here will come from the mix- 
ing of the afterglow contribution with the brighter prompt 
emission component. This mixing will also probably pre- 
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Fig. 4. Early afterglows for a variable source (£em(£) being given by Eq.(13)). Left panel: afterglow computed with 
on- axis emission only compared to the case of a constant source (dotted line). Right panel: afterglow computed with 
off-axis emission included (full line) and geometrical delay At = R/2cT 2 (dashed line). A uniform external medium of 
density n = 1 cm~ 3 and a redshift z — 2 have been assumed. 



vent an unambiguous detection of the imprint of source 
variability on the X-ray afterglow (Fig. 4). In the visible, 
where the burst prompt emission is weak and probably 
negligible (see however the recent RAPTOR observations 
of GRB 041219a and GRB 050820a (Vestrand et al., 2005, 
2006)), the EMM predicts a brighter afterglow for a given 
set of parameters e e , eg, n or A*. But in real afterglows 
these parameters are not known a priori and deciding be- 
tween models will be tricky. Polarization properties of the 
burst prompt emission (Lyutikov, 2004) when they be- 
come more easily accessible may provide clearer evidence. 

A last interesting point concerning the EMM is related 
to the shallow part observed by SWIFT in many X-ray 
afterglows. The light curves in Fig. 3 show that the EMM 
indeed predicts an initially flat region in the early X-ray 
afterglow. However this flat region does not last more than 
the period of source activity (150 s in observer time in 
Fig. 3). It would then extend to 10 4 s (or more) only if the 
source can remain active for that duration, as was also 
suggested for the standard model (Zhang et al., 2005). 



Vestrand, W.T., Wozniak, P.R., Wren, J.A. et al., 2005, Nature 
435, 178. 

Vestrand, W.T., Wren, J.A., Wozniak, PR. et al., 2006, to 

appear in Nature. 
Woods, E., & Loeb, A. 1999, ApJ, 523, 187. 
Zhang, B., Fan, Y.Z., Dyks, J. et al. 2006, ApJ 642, 354. 



References 

Blandford, R.D. & McKee, C.F. 1976, PhFl 19, 1130. 
Daigne, F., & Mochkovitch, R. 1998, MNRAS, 296, 275. 
Granot, J., Piran, T. & Sari, R. 1999, ApJ, 513, 679. 
Jakobsson, P., Levan, A., Fynbo, J.P.U. et al., 2006, A&A 447, 
897. 

Lyutikov, M. & Blandford, R.D., 2003, |astro-ph/0312347l 
Lyutikov, M., 2004, 35th COSPAR Meeting, |astro -ph/0409489 
Sari, R., Piran, T. & Narayan, R., 1998, ApJ 497, L17. 



